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ABSTRACT 

The X-ray luminosity-temperature relation for nearby T ~ 3.5- 10 keV clusters is rederived using new ASCA 
temperatures (Markevitch et al. 1998a) and ROSAT luminosities. Both quantities are derived by directly excluding 
the cooling flow regions. This correction results in a greatly reduced scatter in the Lx~T relation; cooling flow 
clusters are similar to others outside the small cooling flow regions. For a fit of the form Lbd oc T a , we obtain 
a = 2.64 ± 0.27 (90%) and a residual rms scatter in logLbd of 0. 10. The derived relation can be directly compared 
to theoretical predictions that do not include radiative cooling. It also provides an accurate reference point for 
future evolution searches and comparison to cooler clusters. 

The new temperatures and Lx~T relation together with a newly selected cluster sample are used to update the 
temperature function at z ~ 0.05. The resulting function is generally higher and flatter than, although within the 
errors of, the previous estimates by Edge et al. (1990) and Henry & Arnaud (1991, as rederived by Eke et al. 1996). 
For a qualitative estimate of constraints that the new data place on the density fluctuation spectrum, we apply the 
Press-Schechter formalism for J7o = 1 and 0.3. Forfio= 1, assuming cluster isothermality, the temperature function 
implies erg = 0.55 ±0.03, while taking into account the observed cluster temperature profiles, a% = 0.51 ±0.03, 
consistent with the previously derived range. The dependence of cr% on ilo differs from earlier findings, because 
of our treatment of the slope of the fluctuation spectrum n as a free parameter. For the considered values of f2o, 
n = -(2.0-2.3) ± 0.3, somewhat steeper than derived from the earlier temperature function data, in agreement 
with the local slope of the galaxy fluctuation spectrum from the APM survey (Baugh & Gaztanaga 1996), and 
significantly steeper than the standard CDM prediction. 

Subject headings: Cooling flows — cosmology — galaxies: clusters: individual — intergalactic medium — 
X-rays: galaxies 



1. INTRODUCTION 

Clusters exhibit a correlation between their X-ray luminosity 
and temperature, approximately Lboi oc T 3 (Mushotzky 1984; 
Edge & Stewart 1991; David et al. 1993; Fukazawa 1997). 
However, models of cluster evolution that assume no segrega- 
tion of dark and gaseous matter predict Lboi oc T 2 (e.g., Kaiser 
1986; Navarro, Frenk, & White 1995). This discrepancy sug- 
gests that the gas is distributed and evolves differently from the 
dark matter. One of the proposed reasons for this difference is 
preheating of the intracluster gas by an agent other than gravity, 
such as supernovae, at some high redshift (Kaiser 1991; Evrard 
& Henry 1991). Significant energy injection to the gas should 
have been made at the time of its enrichment with heavy el- 
ements (e.g., David, Forman, & Jones 1991; Loewenstein & 
Mushotzky 1996), although neither the amount nor the epoch 
of such injection is certain. Models that allow preheating suc- 
ceed in approximately reproducing the observed Lx~T slope 
(e.g., Navarro et al. 1995). If the gas was preheated and later 
subjected to merger shocks, cooler clusters and groups are ex- 
pected to have a steeper Lx~T dependence than hot clusters 
(e.g., Cavaliere et al. 1997), for which there is some observa- 
tional indication (Ponman et al. 1996). Also, preheating should 
make the redshift evolution of the Lx~T relation slower or 
absent (Evrard & Henry 1991). Indeed, observations exclude 
any dramatic evolution out to z — 0.3-0.5 (Tsuru et al. 1996; 
Mushotzky & Scharf 1997), although this may also be expected 
in low £1() models without preheating (e.g., White & Rees 1978; 
Eke, Navarro, & Frenk 1998). 



The previously reported Lx~T relation at both low and high 
redshifts exhibits a large scatter which precludes accurate deter- 
mination of its exact shape and evolution. Fabian et al. (1994) 
noted that the scatter is mostly due to the strong cooling flow 
clusters. Cooling flow regions, present in more than a half of all 
low-redshift clusters (e.g., Edge, Stewart, & Fabian 1992), are 
governed by a different physics than the bulk of the cluster gas, 
while in extreme cases emitting most of the X-ray luminosity 
and strongly biasing the wide-beam temperature measurements 
(e.g., Allen 1998). Cluster evolution models mentioned above 
predict global properties of the cluster gas and generally ignore 
the presence of cooling flows, thus it is not obvious that their 
predictions can be meaningfully compared to the observations. 
In principle, radiative cooling can be included in the models, 
but it appears that with greater certainty its effects can be sepa- 
rated in the data, which we attempt in this paper. 

The Lx~T relation is a necessary tool for derivation of the 
cluster temperature function using a flux-limited sample. For 
low-redshift clusters, this function was previously obtained by 
Edge et al. (1990, hereafter E90) and Henry & Arnaud (1991, 
hereafter HA; see correction in Eke, Cole, & Frenk 1996). 
Since the gas temperature is linked to the total cluster mass, 
the temperature function provides information on the spectrum 
of the cosmological density fluctuations (e.g., HA), as well as 
on with additional data from higher z (e.g., Henry 1997 and 
references therein). The least certain link in this line of argu- 
ment is the conversion from the observed temperature to the 
cluster mass, which has usually been made assuming cluster 
isothermality or near-isothermality. However, recent ASCA re- 



Also Space Research Institute, Russian Academy of Sciences 



1 



2 



MARKEVITCH 



suits (e.g., Markevitch et al. 1998a, hereafter M98) indicate 
considerable spatial temperature variations which affect the 
temperature-mass relation. In addition, the temperature func- 
tion itself is affected by cooling flows via sample selection, in- 
dividual temperature errors, and the L x -T relation. All these 
effects should be corrected for an accurate comparison with the- 
oretical models. 

In this paper, we use the new ASCA cluster mean tempera- 
tures obtained by directly excluding cooling flow regions (M98) 
and similarly corrected luminosities for a cluster sample se- 
lected from the new ROSAT All-Sky dataset. With these data, 
we rederive the low-redshift L x -T relation and the temperature 
function which can be directly compared to theoretical predic- 
tions as well as to the oncoming high quality data from higher 
redshifts. An accurate comparison with models is best per- 
formed by detailed simulations and is left for future work; for 
a qualitative estimate of the constraints that the new data place 
on the density fluctuation spectrum, we apply a Press-Schechter 
formalism for Qq = I and 0.3. We use Ho = 100 h km s" 1 Mpc" 1 ; 
confidence intervals are one-parameter 90% unless stated oth- 
erwise. 

2. THE SAMPLE 

The sample is selected from the ROSAT All-Sky Survey 
(RASS) Abell cluster list (Ebeling et al. 1996) plus three 
known bright non-Abell clusters, with 0.04 < z < 0.09 and 
/(0.1-2.4keV) > 2 x 10" 11 erg s^cm" 2 . Because the RASS 
fluxes have some uncertainty, pointed ROSAT PSPC or Einstein 
IPC observations of all clusters with the RASS fluxes greater 
than 3/4 of the above limit were analyzed to derive accurate 
fluxes and luminosities. The rederived fluxes correspond to the 
aperture of r = 1 h' 1 Mpc (which includes practically all the 
flux) centered on the X-ray emission centroid, and exclude all 
contaminating sources. Below, these fluxes and corresponding 
luminosities are referred to as "total". 

Cooling flows were excised from fluxes and luminosities al- 
ready at the sample selection step, for consistency with the 
analysis below. A 50 h~ l kpc radius contains most of the cool- 
ing flow emission in nonextreme cooling flow clusters such as 
those in our sample. Therefore, to do the excision in a uni- 
form manner, for all clusters, regions of 50 h~ l kpc radius cen- 
tered on the main brightness peak were masked, and the result- 
ing fluxes and luminosities were multiplied by 1 .06 to account 
for the flux inside the masked region assuming an average (3- 
model for the cluster X-ray brightness, Sx oc (1 + r 2 /a 2 : )" 3/3+1/ ' 2 
with a x = 125/i _1 kpc and j3 = 0.6 (Jones & Forman 1984). 
We refer to these fluxes and luminosities as "corrected". For 
clusters with strong cooling flows (e.g., A478, A780, A2597), 
this correction significantly reduces the flux; for clusters lack- 
ing cooling flows, the flux is (by design) practically unchanged. 
Note that this correction procedure is only approximate and is 
not indended to result in the exact removal of the cooling flow 
contribution, which obviously would require a detailed model- 
ing. However, its simplicity and model-independence enables a 
straightforward comparison to the current simulations that lack 
resolution to reproduce cluster high-density regions. 

Use of more accurate fluxes from pointed observations has 
resulted in adding A3391 and removing A3532 from the origi- 
nal sample, while excision of cooling flows has resulted in re- 
moving A133 and A2597 (with the latter still being used for 
the L x -T relation, as it has been already analyzed). Because 
cooling flows only increase the total flux, with high certainty no 
clusters from the Ebeling et al. (1996) list that satisfy our crite- 



ria were missed. Table 1 lists the final full cluster sample with 
total and corrected luminosities and corrected fluxes. Accuracy 
of the measured fluxes and luminosities is around 5% due to 
the uncertainties in ROSAT calibration, background, Galactic 
absorption, etc. 

TABLE 1 
Full Cluster Sample 



Name 


z 


T 6t 


J corr 




J c 
L x 


Hoi 


fx 


Aoj 


U.UjZ 


A i _i_n o 
O.l ±U.Z 


o.y ±u.4 


7 TO 
Z.ZZ 


1 ft ft 

1.9U 


A ZA 

4.J4 


A C/1 
0.j4 


a i in 

Aiiy 


ft (\aa 
U.U44 


c q i n a 
J.5 ±U.O 


j.o ±u.j 


ft 

U.OJ 


ft CQ 


1 QC 
1 .OJ 


AAA 
4. 14 


A399 


0.072 


7.4 ±0.7 


7.0 ±0.4 


1.66 


1.71 


4. 12 


3.07 


a /im 


U.U / 4 


Q 7 Pi C 
O.J ±U.J 


q n a /i 
o.U ±U.4 


7 Oft 

z.yu 


O Q7 

Z.o / 


7 ACi 

1 Ay 


/I Q7 

4.yj 


A /1TO 

A4/0 


U.Uoo 


*7 i _i_n a 
1 A ±U.4 


Q A +0.8 
8 ' 4 -1.4 


4.1 j 


1 7 1 

j.Zl 


O A7 

o.oZ 


7 on 
3.5 / 


A A /I A * 

Ao44 T 


U.U / 1 


"7 1 A A 

/.I ±U.o 


"7 ft J_ft O 

/.y ±u.o 


7 ftT 
Z.U/ 


1.0/ 


A OQ 
4.0J 


7/1/1 

j.44 


A / J4 


U.Uj4 


q n i o c 
y.u ztu.j 


a c, +0.7 
y " 5 -0.4 


1 Q 

z.iy 


Z.ZZ 


0.44 


/.u / 


A TQft 


U.U j / 


Q Q A 7 

3.5 ±U.Z 


4. j ±U.4 


1 AO 

l.oo 


1 1 "7 
1.1 / 


7 17 

Z. 1Z 


a /i ft 
3.40 


A 1644 


ft C\A Q 

U.U40 


A 1 I A q A 

4. / ±U.o a 


5.0 e 


A ft/i 
u.y4 


0.93 




3.79 


a i Acft 
AIojO 


n ftoc 
U.Uoj 


C A I A 
J.O ±U.O 


A T A O 

o. / ±U.o 


1.50 


l.OJ 


A 11 


2.39 


A 1 AC 1 

AIoj 1 


ft ftoc 
U.Uoj 


A 7 A C 
O.J ±U.J 


A 1 I A ,1 

0.1 ±U.4 


1 ftl 

1.9/ 


1 Oft 

l.oU 


A ftT 
4.0/ 


7 ac 
Z.jj 


A 1 T2A 
Al /JO 


U.U40 


1 C A /l 

J.J ±U.4 


J.J ±U.4 


ft eft 
U.jU 


ft Cft 

U.jU 


ft 07 

U.OJ 


7 1 Q 

Z.lo 


A 1 TOC 

Al /9j 


U.Uoz 


a n i a 7 
O.U ±0.3 


1 Q -L- 1 ft 

/.o ±1.U 


7 Aft 

Z.oU 


1 ftO 

1.9o 


C ftft 

j.uy 


A "70 

4. /o 


a 7ft70 
AZUZy 


U.U / / 


Q "7 A T 
5. / ±U.J 


o 1 i i a 
y.i ±i.u 


4.Z0 


J.JO 


ft A ft 

y.4u 


C 7ft 

j.jy 


AZUOj 


n mi 
U.U / z 


c A I n a 
J. 4 ±U.j 


C C A A 

J.J ±U.4 


1 1Q 

i.zy 


1 07 
l.Z / 


7 AA 
Z.00 


O O/l 
Z.Z4 


A2142 


0.089 


8.8 ±0.6 


a 7 +1.5 

9 - 7 -l.l 


5.22 


4.83 


14. 19 


C AO 
J.OO 


AZZjO 


U.U J 5 


"7 < A /I 

/.J ±U.4 


0.0 ±U.4 


Z.l J 


1 1Q 

Z. lo 


C ftA 
J.UO 


A ft1 

o.ui 


AZJiy^ 


U.U jo 


ft 1 J_ft "7 

y.z ±u. / 


O O A c 

0.0 ±U.J 


7 CH 
J.O/ 


7 CA 
J.04 


1 Ct f^A 

lU.Uo 


1 1 ftft 
i i.uy 


A O C CO 

AZjoy 


ft ft/17 

U.U4Z 


q 7 +2.2 j 

3.7 _ 12 d 


j.j e 


ft <1 

U.Jl 


ft AA 
U.44 




7 1Q 
Z.JO 


AZjy /j 


n noc 
U.Uoj 


i.o ±0.2 


/I /I +0.4 


1 7 -a 
1. /j 


ft ftO 

U.98 


1 O 1 

l.ol 


1 7 A 
l.ZO 


A 7AC7 
AZOj / 


U.U4U 


i "7 i 


3. / ±U.3 


ft /17 

U.4 / 


ft AO 
U.4Z 


ft 71 

U. / 1 


Z.44 


A3112 


0.070 


4.7 ±0.4 


c o +0 7 
J - J -1.0 


1.64 


1.15 


2.34 


2.16 


A3158 


0.059 


5.5 ±0.6 d 


6.0 e 


1.41 


1.38 




3.68 


A3266 


0.055 


7.7 ±0.8 


8.0+0.5 


1.82 


1.83 


4.78 


5.76 


A3376 


0.046 


4.3 ±0.6 


4.0 ±0.4 


0.56 


0.57 


0.99 


2.44 


A3391 


0.054 


5.7 ±0.7 


5.4 ±0.6 


0.74 


0.74 


1.52 


2.39 


A3395 


0.050 


4.8 ±0.4 


5.0 ±0.3 


0.79 


0.81 


1.60 


3.05 


A3558 


0.048 


5.5 ±0.3 


5.5 ±0.4 


1.74 


1.68 


3.51 


6.89 


A3562 


0.050 


3.8 ±0.9 d 


4.8 e 


0.90 


0.87 




3.21 


A3571 


0.040 


6.9 ±0.3 


6.9 ±0.2 


1.97 


1.83 


4.36 


10.78 


A3667 


0.053 


7.0 ±0.6 


7.0 ±0.6 


2.18 


2.22 


5.34 


7.36 


A4059 


0.048 


4.1 ±0.3 


4.4 ±0.3 


0.81 


0.66 


1.22 


2.67 


Cyg A* 


0.057 


6.5 ±0.6 


6.1 ±0.4 


2.70 


2.01 


4.46 


5.79 


MKW3S 


0.045 


3.5 ±0.2 


3.7 ±0.2 


0.75 


0.54 


0.91 


2.47 


Triang* 


0.05 1 


9.5 ±0.7 


9.6 ±0.6 


3.12 


3.13 


9.13 


11.21 



NOTES : Temperatures are in keV; luminosities are in units of 1 44 /T 2 erg s _1 ; 
fluxes are corrected for absorption and in units of 10~" erg s~' cm -2 . Subscript 
X denotes 0.1-2.4 keV energy band; bol is for bolometric. Clusters with * are 
at < 20° and the flux of the cluster with f is below the limit; these clusters 
are used for the Lx — T relation only. 

a Single-temperature fits and cooling flow-corrected emission-weighted tem- 
peratures with 90% errors are from M98, except those marked d which are from 
David et al. (1993) and e which are estimated from the Lx — T relation derived 
here. 

b Total luminosities within r < 1 hr l Mpc (including cooling flow regions). 
c Corrected quantities (cooling flow regions excised). 

For all clusters in the sample except 4, ASCA wide-beam as 
well as emission-weighted cooling flow-excluded temperatures 
are derived in M98 and references therein. These temperatures 
correspond to averaging over r = 0.6- 1 h~ x Mpc apertures, de- 
pendent on the individual observation details; because the total 
cluster emission is dominated by bright central regions, such a 
difference in radial coverage is unimportant. For the remaining 
4 clusters without the ASCA data, temperatures are estimated 
from the L x -T relation derived in §. Non-ASCA wide-beam 
temperatures for these clusters exist (David et al. 1993) and 
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FIG. 1. — Luminosities in the 0.1-2.4 keV band vs. temperatures from Table 1. Panel (a) shows total luminosities within r = Iff 1 Mpc and 
single-temperature fits to the overall ASCA spectra, while panel (b) shows corrected luminosities with cooling flow regions excised (see text) 
vs. emission-weighted temperatures derived from the temperature maps excluding cooling flow components. Triangles mark clusters with strong 
cooling flows, circles are clusters with weak or no cooling flows. Dashed lines are respective best-fit power laws from Table 2. 



are consistent, within their relatively large uncertainties, with 
the estimated values. We will use our estimates for the tem- 
perature function, since those measured temperatures are not 
cooling flow-corrected. For clusters with ASCA temperatures, 
Table 1 also lists cooling flow-excised bolometric luminosities, 
estimated from the corrected 0. 1-2.4 keV luminosities using the 
corrected temperature and assuming isothermality. Our redshift 
interval, imposed by the ASCA analysis constraints (see M98), 
effectively excludes clusters with T <; 3.5 keV. 

For the temperature function determination, as opposed to 
the Lx~T relation, completeness of the sample is important. 
Therefore, as in previous analyses, the \b\ < 20° area is ex- 
cluded from the temperature function subsample. Ebeling et al. 
(1996) find that their originally optically selected list is already 
noticeably incomplete in the interval 20° < \b\ < 30°. Exclud- 
ing this area does not change our results, so this small effect is 
ignored. Note that our sample, which consists of bright, mod- 
erately nearby clusters with a flux limit a factor of 4 above the 
completeness limit of the Ebeling et al. (1996) list, is less likely 
to be affected by the above effect than their full list. In §, we 
confirm the completeness of our sample using the luminosity 
function from the purely X-ray selected RASS cluster catalog 
(Ebeling et al. 1997; the catalog itself is unavailable at the time 
of this writing). 

The subsample of 30 clusters used for the temperature func- 
tion has a median temperature and redshift of about 6 keV and 
0.054, respectively. It overlaps to a large degree the samples 
of E90 and HA, both selected by 2-10 keV flux, but differs 
from them mainly at low temperatures and fluxes. Unlike our 
sample, the E90 dataset is not limited by redshift, which im- 
proves statistics by including several hot, distant clusters, but 
has the disadvantage of complicating the evolutionary studies. 
Our sample has a lower effective flux limit than HA and there- 
fore has more clusters in the corresponding temperature range. 
The main improvement upon previous datasets is, of course, 
the higher-quality imaging and temperature data, including spa- 
tially resolved temperature measurements, available to us. 



3. THE Lx-T RELATION 

Table 1 shows that for strong cooling flow clusters, the cool- 
ing flow-corrected temperatures are significantly higher than 
the wide-beam values, and that the small central regions con- 
tain a significant fraction of the total X-ray luminosity. Below 
the luminosity-temperature relation in its traditional meaning, 
using wide-beam temperatures and total luminosities, is com- 
pared to the one that corresponds only to the main cluster gas, 
excluding cooling flows in clusters that have ones. 
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FIG. 2. — Same as Fig. lb (luminosities vs. temperatures both cor- 
rected for cooling flows) but for bolometric luminosities. 

Figure 1 shows luminosities in the 0. 1-2.4 keV band vs. tem- 
peratures for clusters with ASCA data, including and excluding 
the cooling flow regions. Corrected bolometric luminosities 
vs. corrected temperatures are plotted in Fig. 2. These L-T 
relations are fit by power laws of the form L = A 6 r 6 ", where 
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T(, = T/6 ke V, using the bisector modification of the Akritas & 
Bershady (1996 and references therein) linear regression algo- 
rithm that allows for intrinsic scatter and nonuniform measure- 
ment errors and treats L and T symmetrically. The confidence 
intervals are evaluated by bootstrap resampling (e.g., Press et 
al. 1992) while simultaneously adding random measurement er- 
rors to the temperatures and 5% errors to the luminosities (the 
results do not depend on the latter). Table 2 lists the best-fit 
parameters together with the residual rms scatter in logL for a 
given T, <j\ og i (the scatter along the T axis is <j\ og T — ciogz,/oO. 
We also include a fit for the cooling flow-corrected tempera- 
tures vs. total luminosities, which may be useful, for example, 
to relate theoretical predictions to the fluxes of poorly resolved 
distant clusters. Note that the 0.1-2.4 keV gas emissivity is 
only weakly dependent on temperature in the considered tem- 
perature range, so the slope a of the L^i - T relation should be 
greater than that for the 0. 1-2.4 keV band by about 0.5 without 
any corrections. 



TABLE 2 
Fits to L-T Relations 



Band 


Data 


A 6 


a 


flogi 


keV 




10 44 r 2 ergs-' 






0.1-2.4 


Uncorr. T, total L 


1.71 ±0.21 


2.02 ±0.40 


0.181 


0.1-2.4 


Corr. T, total L 


1.57 ±0.17 


2.09 ±0.29 


0.133 


0.1-2.4 


Corr. T, corr. L 


1.41 ±0.12 


2. 10 ±0.24 


0.104 


Bol. 


Corr. T, corr. L 


3.11 ±0.27 


2.64 ±0.27 


0.103 



Figure 1 and Table 2 show that the scatter in the L x -T rela- 
tion is greatly reduced when cooling flows are excised and the 
quantities representing only the main cluster gas are compared. 
An example of the two largest deviations around 3.7 keV in 
Fig. la, which are A780 (Hydra A) and A2597, shows that the 
improvement is due in equal degrees to the corrections of the 
temperature and luminosity. The same can be seen in the last 
column of Table 2. Such a strong effect of the cooling flows on 
the wide-beam temperatures was not anticipated before ASCA, 
and the cooling flow luminosity alone is insufficient to account 
for the scatter. This led Fabian et al. (1994) to suggest that 
cooling flow clusters are different from others in their global 
gas properties. Our results show that most of the scatter in the 
L x -T relation is due to localized, r <^ 50/T 1 kpc regions with 
cooling flows. On the other hand, after the exclusion of these 
regions, some residual scatter remains, which indeed is likely to 
be due to the different gas distributions — after the corrections, 
the largest deviation at 9.5 keV in Fig. lb is a strong merger 
A754. Note that a small fraction of the scatter is due to the 
temperature measurement errors. 

The derived L x -T relation and its intrinsic scatter can 
be compared to theoretical and numerical predictions that at 
present do not model radiative cooling adequately. A self- 
similar model of cluster growth (Kaiser 1986), supported by 
simulations that include only gravity (e.g., Navarro et al. 1995), 
predicts Lboi °c T 2 . A disagreement of this prediction and the 
observed steeper slope has long been noted; our new Lboi _ T 
relation is also significantly steeper. To explain this disagree- 
ment, it was proposed that hypothetical preheating of the intra- 
cluster gas, e.g., by supernovae, effectively reduces the central 
gas densities with a stronger effect on the cooler clusters (e.g., 
Kaiser 1991; Evrard & Henry 1991; Navarro et al. 1995; Cava- 
liere et al. 1997). Assuming that, after initial preheating, gas in 
the cluster cores remains on the same adiabat as clusters evolve, 
Evrard & Henry predict a scale-free relation L bo i oc T 2JS . How- 



ever, cluster merging should result in significant increase of the 
gas specific entropy due to shock-heating, as is indeed observed 
(e.g., Markevitch et al. 1998b). Cavaliere et al. allow the gas in- 
side the cores to be shock-heated and mixed in addition to an 
early injection of the energy equivalent to 0.5-0.8 keV per gas 
particle. Their predicted L DO i - T slope changes from 5 on the 
group scale to 2 on the hottest cluster scale, being ~ 2.3 in the 
temperature interval covered by our data. They also predict a 
scatter in the Lboi - T relation of about cri ogi ~ 0.13 due to the 
differing merging histories of individual clusters. Considering 
the number of adjustable quantities in such modeling, these pre- 
dictions are quite close to our observations. Although our data 
are adequately described by a single power law, with such a 
small scatter it is quite feasible to measure the predicted steep- 
ening of the L-T slope toward lower temperatures, if compa- 
rably accurate data at those temperatures are obtained. Such a 
measurement would provide interesting information on the en- 
ergetics of the intracluster gas. 

The observed tight L x -T correlation has other implications. 
M. Arnaud (1997, private communication) used the small scat- 
ter of L x for clusters lacking cooling flows to constrain the 
cluster-to-cluster baryon fraction variations. We also note for 
illustration that our scatter in low-redshift L x is equivalent to 
0.25" 1 , which is comparable to ~ 0.19 m for Type la super- 
novae used for cosmological distance estimates (Perlmutter et 
al. 1997). However, clusters are not expected to remain stan- 
dard candles at all redshifts. 

After this paper has been submitted for publication, we re- 
ceived a preprint by Allen & Fabian (1998) who employed a 
different approach to evaluate the effect of cooling flows on 
the L x -T relation for hot clusters. They modeled the over- 
all cluster spectra by a sum of the isothermal and cooling flow 
components and took the temperature of the isothermal com- 
ponent as the true average cluster temperature. This compares 
to the M98 approach of using spatially resolved spectra to lo- 
calize cooling flow regions, obtain temperature maps outside 
these regions, and use these maps to calculate average temper- 
atures. Allen & Fabian find generally similar temperature bias 
for those clusters studied in both works, although the exact tem- 
perature values are sometimes different due to the difference in 
approach. Unlike our result, the Allen & Fabian's corrected 
L x -T relation still has considerable scatter and shows differ- 
ence between clusters with and without cooling flows. This is 
because their luminosities include cooling flow regions, while 
we excise these regions when calculating the corrected lumi- 
nosities. Nevertheless, the slope of the corrected Lboi oc T a 
relation obtained by Allen & Fabian, a = 2.3^ 3, is in good 
agreement with our result. 

4. THE TEMPERATURE FUNCTION 

4.1. Derivation 

The temperature function (TF), which is the number of clus- 
ters per unit comoving volume and unit temperature interval, 
dN/dT, or the number above a given temperature, N(> T), 
can be derived by co-adding all clusters with weights 1 /V(T). 
Here V(T) is the maximum comoving volume within which 
a cluster with a given temperature could have been detected 
above the flux limit of our sample, and T is a cooling flow- 
corrected temperature. Since the sample is flux-limited rather 
than temperature-limited, an L x - T relation has to be used to 
calculate the volume as a function of temperature. The sample 
was selected using the fluxes with excised cooling flows and we 
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can take advantage of the tight Lx~T relation obtained above. 
The volume as a function of temperature was calculated by in- 
tegrating, over the redshift interval 0.04-0.09, volume elements 
with weights equal to the probability of a cluster with a given 
temperature to have a detectable flux, assuming a log-Gaussian 
scatter in the Lx~T relation given in Table 2. For a purely flux- 
limited sample, a symmetric Lx~T scatter such as that results 
in a uniform increase of the volume by a factor of about 1 .06, 
compared to no scatter.^] In the presence of the redshift con- 
straints, the scatter results in the volume correction between a 
factor of 1.26 at our minimum T = 3.5 keV and 0.88 at 6.5 keV. 
The sample is essentially volume-limited above ~ 7 keV. 

Note that this approach to the derivation of the TF differs 
from the one employed, e.g., in HA and other works, in which 
individual cluster luminosities (rather than temperatures) are 
used to determine the volume within which a given cluster 
would have a detectable flux. The latter method is entirely ap- 
propriate for the derivation of a luminosity function, but it is 
not optimal for the TF, because clusters with a given tempera- 
ture have a wide range of luminosities. This method does result 
in an unbiased, on average, TF estimate in the limit of the large 
number of clusters. In reality, however, the number of clusters 
is always limited, and the large intrinsic scatter in the Lx~T 
relation causes bias in the median and most probable TF es- 
timates, especially at the highest (lowest) temperatures where 
the TF is estimated using only one or two hottest (coolest) clus- 
ters in the sample. To quantify this bias for a purely flux-limited 
sample, we performed a simple simulation, scattering randomly 
in Euclidean space clusters with luminosities distributed with 
an rms Gaussian scatter in \ogLx of 0.28 (which corresponds 
to the Lx~T scatter in HA). For those simulated clusters above 
the flux limit, search volumes were calculated using individual 
luminosities, V oc L 3 / 2 . The resulting distribution of the calcu- 
lated weights 1 /V is strongly skewed; the median derived TF 
is only 0.65 of the average (true) value at a temperature of the 
single hottest cluster in the sample, while the most probable TF 
value is only 0.32 of the average. The reason is that for any 
temperature, there are many more overluminous clusters in a 
flux-limited sample than underluminous ones (for the adopted 
value of the Lx~T scatter, the ratio is 5:1).[] The TF deriva- 
tion method employed in this paper is free of this small-number 
bias. Its possible disadvantage is the reliance on the validity of 
the Lx~T relation at all temperatures, including the extreme 
ones where the number of clusters in the Lx~T sample is small. 
However, this does not apply to our sample, since it is volume- 
limited above ~ 7 keV, which means the search volume for hot 
clusters is simply the total volume between 0.04 < z < 0.09. 

The uncertainty of the resulting temperature function is eval- 
uated by bootstrap resampling, assuming a Poissonian distribu- 
tion of the total number of clusters in the sample and adding 
random Gaussian measurement errors to the individual temper- 
atures at each resampling. The latter is a minor contribution to 
the overall uncertainty dominated by that of the small number 
of clusters. 

The resulting cumulative temperature function is plotted in 
Fig. 3, overlaid in panel (a) on the previously derived functions 
of E90 and HA (the latter is taken from Eke et al. 1996). The 
figure shows that the overall effect of our corrections for the 
presence of cooling flows at each stage of the derivation is not 



large — there is agreement with the earlier E90 and HA func- 
tions within their ler uncertainties (not shown). The cooling 
flow correction on average results in a small shift of the TF to- 
ward higher temperatures. An apparent (not highly significant) 
disagreement between our and HA results at high temperatures 
is mostly due to (a) the fact that our sample is volume-limited 
while the HA sample is flux-limited, (b) a possible incomplete- 
ness of the HA sample (e.g., it does not include the hot cluster 
A2163 which appears to satisfy the selection criteria), and (c) 
the possibility of the small-number bias described above. 

To assess the completeness of our sample, particularly at low 
temperatures where the Abell catalog may miss objects even at 
high Galactic latitudes (e.g., MKW3S), we use the luminosity 
function obtained from the purely X-ray selected RASS cluster 
catalog (Ebeling et al. 1997) and convert it into the temperature 
function using the Lx~T relation and its scatter. The relevant 
Lx~T relation was derived using the total, RASS -measured 
(Ebeling et al. 1996) luminosities and corrected temperatures 
of our clusters. The resulting relation is similar to line 2 in Ta- 
ble 2 but has a larger scatter due to the RASS flux uncertainties. 
The resulting temperature function, shown in Fig. 3a, is in ex- 
cellent agreement with our result, indicating that no significant 
number of cool clusters is missing. The luminosity function 
was derived from a much larger cluster sample than ours so its 
uncertainty is negligible. 

To obtain an analytic fit to the differential temperature func- 
tion, we use fine temperature binning and maximum likelihood 
minimization assuming a Poissonian distribution of the cluster 
number in each bin (bins are allowed to have zero clusters; the 
results are independent of the bin size). The errors are calcu- 
lated by the bootstrap algorithm outlined above. For a power 
law of the form dN/dT = A 6 T^ a , we obtain, assuming f^o = 1, 
A 6 = 1.90 ±0.54 x 10- 7 /i 3 Mpc- 3 keV _1 and a = 4.2 ± 0.7. A 
slightly better fit is obtained for an exponential Aeexp[-(r- 
6 keV)/7;] with A 6 = 2.38 ±0.68 x 10- 7 /i 3 Mpc" 3 keV _1 and 
r„ = 1 .6 ± 0.4 keV. These models are shown in cumulative form 
in Fig. 3b. 

4.2. Constraining the Fluctuation Spectrum 

The cluster temperature function was used to obtain con- 
straints on the cosmological density fluctuation spectrum by 
many authors (the incomplete listing includes HA; Kaiser 1991; 
White, Efstathiou, & Frenk 1993; Viana & Liddle 1996; Eke 
et al. 1996; Colafrancesco, Mazzotta, & Vittorio 1997; Pen 
1998) applying the Press & Schechter (1974) formalism sup- 
plemented by simulations. Our temperature function is slightly 
higher and flatter than the previous estimates, and its accuracy 
is better. In addition, Markevitch & Vikhlinin (1997) and M98 
recently showed that the observed cluster nonisothermality con- 
siderably reduces the X-ray mass estimate for a given emission- 
weighted temperature, compared to the previously made esti- 
mates. An accurate comparison of cosmological models to the 
data requires hydrodynamic simulations to take into account 
the nonequilibrium state of many clusters. However, we can 
perform a qualitative estimate of the effect of the above ob- 
servational updates on the fluctuation spectrum constraints by 
applying the Press-Schechter formalism for Oo = 1 an d = 0.3 
(A = 0), and comparing the results to similar previous analyses. 

2 The reason is geometric; the redshift shell [z, z + Az] containing the overluminous clusters has greater volume than the shell [z— Az, z] containing the underlumi- 
nous ones. 

3 In principle, this may also affect the validity of the Lx-T relation derived from a flux-limited sample. However, for our sample, selected by the corrected fluxes 
and with a much smaller scatter in the corrected Lx~T relation, this effect is small. 
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FIG. 3. — Cumulative temperature function. Panel (a) shows our result (assuming J7o = 1) with its 68% error band overlaid on the earlier deriva- 
tions of HA and E90, as well as a function obtained from the luminosity function for a purely X-ray selected sample (Ebeling et al. 1997) using the 
Lx-T relation (see text). Panel (b) shows fits to our data: Press-Schechter fits for n = -1 (standard CDM) and n = -2.3 (best-fit for Qo = 1), and 
empirical power law and exponential fits (see text). 



A full description of the formalism can be found in the refer- 
ences above, so only the relevant definitions are given here. 

We assume that the power spectrum of linear density fluctua- 
tions on the cluster scale is oc k", implying an rms amplitude 
of the linear mass fluctuations in a sphere containing an av- 
erage mass M of a(M) = <r% (M/m 8 )" ( " +3)//6 (l +z) _1 , where m 8 
is the average mass within the comoving r = 8/T 1 Mpc. Our 
clusters are assumed to lie at the median sample redshift of 
0.054, being observed just after their collapse (a reasonable 
assumption for f^o = 1 but not necessarily such for Slo <C 1, 
e.g., White & Rees 1978; nevertheless, we use this assump- 
tion for Hq = 0.3 for comparison with other works). The cluster 
mass within a radius of overdensity A c (taking A c = 178 for 
tt = 1 and A c = 124 for ft = 0.3 as in Eke et al. 1996), is as- 
sumed to scale with the X-ray emission-weighted temperature 
as M = MiofT/lO keV) 3/2 with an rms scatter of 20% as pre- 
dicted by simulations of Evrard, Metzler, & Navarro (1996). 

It is worth mentioning here that the scaling M oc T 3 / 2 follows 
from the same line of argument as the relation Lboi oc T 2 , and 
the latter is not observed. Processes that modify the Lx~T re- 
lation, e.g., preheating, can possibly break the M-T scaling. 
However, X-ray luminosity is very sensitive to the gas density 
distribution, which is more likely to change due to preheating 
than the gas temperature. Indeed, simulations by Metzler & 
Evrard (1998) suggest that the M-T relation holds at least for 
their considered preheating model, while the gas density dis- 
tributions (and therefore the Lx~T relation) are significantly 
modified. The available observational data (masses from grav- 
itational lensing vs. X-ray temperatures) at least do not contra- 
dict the adopted M-T relation (Hjorth, Oukbir, & van Kampen 
1998). X-ray measurements of masses at different temperatures 
(our work in progress) will further test the accuracy of this scal- 
ing. 

For M10 we use two values for each assumed f2o (values in 
parentheses correspond to f2 = 0.3), 13.5 (16.4) x 10 14 /r' Mq 
which assumes cluster isothermality (e.g., Eke et al. 1996) and 
9.1 (9.8) x 10 14 /z _1 Mq which corresponds to the mass pro- 



file of the main cluster of A2256 derived by Markevitch & 
Vikhlinin using the observed temperature profile. Although 
the latter values are based on a single cluster, temperature pro- 
files for many other clusters are similar and the effect on the 
mass is expected to be similar as well (M98). Using these re- 
lations, m% = 6 x 10 14 f2()/i -1 Mq corresponds to T ~ 3-8 keV 
for Slo = 0.3- 1, indicating that the cluster data best constrain 
the density fluctuations at this particular scale. 

The Press-Schechter mass function was converted to the dif- 
ferential temperature function which was then fitted to the data 
as described in §, with a% and n being free parameters. The 
observed TF was derived self-consistently using the two val- 
ues of f^o- The best-fit a% and n are given in Table 3 for the 
two cosmologies and two assumed cluster temperature distribu- 
tions. The best-fit model TF in the cumulative form for f2o = 1 
is shown in Fig. 3b (the two temperature profile cases are al- 
most identical). As was noted in previous works, because the 
abundance of clusters is very sensitive to <j%, the latter quan- 
tity is rather robust to the uncertainty in the number of clusters. 
Therefore, even though our temperature function is about a fac- 
tor of 2 above that of HA used in most previous works, our 
us for the isothermal case and Slo = 1 is only slightly higher 
than the values derived from the corrected HA data (0.50, Eke 
et al. 1996, 1998; 0.53, Pen 1998). Our result for O = 0.3 is 
significantly different from the above two works which report 
<j% — 0.9 for £lo = 0.3. The difference is due to our treatment of 
n as a free parameter; fixing n ~ -1 .3 or -1 .4 as in those works 
and assuming isothermality results in best-fit 08 values close to 
those in the above works, for both values of fio- However, such 
values of n do not fit the data. 



TABLE 3 

Press-Schechter Fits to the Temperature Function 



T PROFILE 


n -- 


= 1 


n a = 


0.3 






n 




n 


Isothermal 


0.55 ±0.03 


-2.3 ±0.3 


0.78 ±0.04 


-2.0 ±0.3 


Observed 


0.51 ±0.03 


-2.3 ±0.3 


0.66 ±0.03 


-2.1 ±0.3 
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FIG. 4. — Circles with la errors show the slope of the APM galaxy 
density fluctuation spectrum (ft oc k") as a function of scale, from 
Baugh & Gaztanaga (1996). Flattening of the observed spectrum at 
k > 0.2 reflects the onset of nonlinear fluctuation growth; dotted line 
is a fit to the de-evolved, linear spectrum. Rectangles shows our slope 
of the linear spectrum (solid line for Qo = 1, dashed line for Qo = 0.3). 
Their vertical sizes are the 90% errors, while their horizontal sizes 
illustrate the approximate k coverage corresponding to our tempera- 
ture range. The rectangles are plotted at k = 1 /r, where r is the radius 
containing the average mass M corresponding to a given temperature. 
A more meaningful average k weighted by the contribution of each 
harmonic to a(M) depends on the exact spectrum, so the rectangle 
positions are rather approximate. 

The derived a% is sensitive to the assumed mass-temperature 
relation. The observed nonisothermality, which modifies this 
relation, acts to reduce the derived value (see Table 3); the scat- 
ter in the relation, included in all our estimates, also reduces the 
us estimate by ~ 0.01. Since some of the observed radial tem- 
perature decline may be due to the incomplete thermalization 
of the gas in the outermost cluster regions (biasing low the X- 
ray-derived total mass), detailed cluster simulations reproduc- 
ing such a decline are necessary for an accurate interpretation 
of the temperature function. In the absence of such simulations, 
the a% values obtained above may be viewed as conservative 
brackets of the true value for the assumed cosmologies. 

The best-fit fluctuation spectrum slope is n = -2.3 and n ~ -2 
for = 1 and fio = 0.3, respectively; the results are very close 
for isothermal and nonisothermal cases (see Table 3). The value 
corresponding to O = 1 is consistent with n~-2 obtained by 
Kaiser (1991) using the E90 data and is marginally steeper than 
n = -(1.7^-35) (68% two-parameter interval) reported by HA 
and n = -1 .8 ± 0.4 by Colafrancesco et al. (1997) using the HA 
data. Our weak dependence of the best-fit n values on fio seems 
to disagree with the Colafrancesco et al. finding; however, it is 
difficult to compare fits to the TFs with rather different observed 
shapes. 

Our temperature interval corresponds to the log range in the 
wavenumber k of only 0.23 (although a comparable interval of 
k contributes to the mass fluctuation a(M) at any M), thus it 
probes a local shape of the fluctuation spectrum. Coinciden- 
tally, the galaxy fluctuation spectrum derived from the APM 
survey, which is the most accurate determination to date (Baugh 
& Efstathiou 1993 and references therein; Baugh & Gaztanaga 



1996), exhibits maximum steepness at this scale with a local 
slope of n ~ -1.9 ±0.3 (68% interval). Figure 4 shows the 
observed APM slope as a function of scale from Baugh & Gaz- 
tanaga together with our measurement. Flattening at k ;> 0.2 
reflects the onset of nonlinear fluctuation growth; the above au- 
thors recover the approximate primordial, linear spectrum (dot- 
ted line in Fig. 4) with n ~ -2 on these scales. The recon- 
structed linear-regime slope as well as the observed slope are 
in agreement with our X-ray results. It is worth recalling that 
the galaxy distribution is linked to the underlying, almost lin- 
ear, density field via an unknown and possibly scale-dependent 
bias. Clusters, on the other hand, represent nonlinearly evolved 
density peaks and their abundance constrains the density spec- 
trum directly, under the assumption of gaussianity of the fluctu- 
ations (and, of course, assuming that the X-ray-measured clus- 
ter masses are not grossly in error). Agreement between two 
such widely different methods is encouraging. Both the optical 
and the X-ray slope are considerably steeper than the standard 
CDM prediction of n ~ -1 at this scale (e.g., Blumenthal et al. 
1984). For illustration of the discrepancy with the X-ray data, 
Fig. 3b shows a temperature function for n = -1. A discussion 
of some promising alternatives (namely, mixed dark matter) can 
be found, e.g., in Baugh & Efstathiou (1993). 



5. SUMMARY AND FUTURE WORK 

Using the new ASCA cluster temperatures and ROSAT lu- 
minosities, both calculated directly excluding the cooling flow 
regions, we obtain the L x -T relation in the 3.5-10 keV temper- 
ature interval with greatly reduced scatter. This result provides 
an accurate reference point (a) for a comparison with nearby, 
cooler clusters and groups that should have a different slope if 
there was preheating, and (b) for evolution studies using the 
oncoming AXAF data which will provide the necessary angular 
resolution for cooling flow excision at higher redshifts. The de- 
rived Lx~T relation can be directly compared to theoretical and 
numerical predictions, most of which do not include a detailed 
treatment of radiative cooling. 

These new data are also used to rederive the nearby clus- 
ter temperature function. This function is interpreted qualita- 
tively by applying the Press-Schechter formalism for fio = 1 
and fio = 0.3, to obtain values of erg assuming the isothermal 
and observed cluster temperature profiles (see Table 3). These 
values should bracket the correct value if some of the observed 
radial temperature decline in clusters is due to the incomplete 
thermalization of the gas in the outer regions. The derived 
slope of the mass fluctuation spectrum at the cluster scale is 
n = -(2.0- 2.3) ±0.3 for considered values of Slo, which is con- 
sistent with the APM galaxy survey at the same scale. Again, 
our temperature function provides a reference for future accu- 
rate evolution studies and for direct comparison with theoretical 
predictions. 

It would be interesting to extend the temperature range of this 
work to lower temperatures, 1-3 keV, where one expects to see, 
for example, the effects of preheating on the Lx~T relation. 
It would also greatly strengthen the constraints from the tem- 
perature function. Such an extension will be possible when X- 
ray selected cluster catalogs become available. Since preheat- 
ing can modify the simple mass-temperature scaling, accurate 
measurements of cluster masses over a range of temperatures, 
including the lowest, will be required for the interpretation of 
such an extended temperature function. 
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